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ABSTRACT 

*n : 

Aims. Hydrogen and helium line spectra are crucial diagnostic features for the quantitative analysis of OB stars. Hybrid non-LTE line- 
CNl ■ formation calculations for these elements have not been discussed thoroughly so far, despite their wide use for analyses of metal line spectra. 
Ji* We compute synthetic spectra based on a hybrid non-LTE approach in order to test the ability of these models to reproduce high-resolution and 
high-S/N spectra of dwarf and giant stars and also to compare them with published grids of non-LTE (OSTAR2002) and LTE (Padova) models. 
' ' Methods. Our approach solves the restricted non-LTE problem based on classical line-blanketed LTE model atmospheres. State-of-the-art 
model atoms and line-broadening theories are employed to model the H and He i/ii spectra over the entire optical range and in the near-IR. 
Results. A comparison with published line-blanketed non-LTE models validates the suitability of the LTE approximation for modelling 
the atmospheric structure of late O to early B-type dwarf and giant stars at metallicities down to (at least) 1/5 x solar. Our hybrid non-LTE 
, synthetic spectra simultaneously match almost all measurable hydrogen and helium lines observed in six test stars over a wide spectral range 
from the Balmer limit to the near-IR, except for only a few well-understood cases. A robust starting point for studies of the metal spectra 
is thus established. Our approach reproduces other published non-LTE calculations, but avoids inconsistencies in the modelling of the He i 
singlets. These have recently been discussed in the literature in the context of O-type stars and we show that they persist in the early B-types. 
Q ■ Our approach improves on published pure LTE models - widely applied for OB star analyses - in many aspects: non-LTE strengthening and 
< the use of improved line-broadening data result in significant differences in the line profiles and equivalent widths of the Balmer and helium 
C/3 lines. Where possible, systematic effects on the stellar parameter determination are quantified, e.g. gravities derived from the Hy wings may 
be overestimated by up to ~0.2 dex for dwarfs at our upper temperature boundary of 35 000 K in LTE. 
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1. Introduction Quantitative analyses of O- and B-type stars can provide 

observational constraints on both stellar and galactochemical 
Massive stars are the main drivers for the dynamical and chem- evolution. In the first case, information on basic stellar param- 
ical evolution of the interstellar medium (ISM), and thus for the eters and mixing with nuclear-processed matter (abundances 
evolution of entire galaxies. They are important sites of nucleo- of the light elements) facilitates empirical evaluation of the 
synthesis and main contributors to the energy and momentum quality of different evolutionary models (see e.g. Maeder & 
budget of the ISM, as sources of ionizing UV radiation, through Mevnet [2000T> . In the second case, information on the spa- 
stellar winds and their final explosion as supernovae. tial distribution of present-day elemental abundances allows 

Galactic metallicity gradients to be derived (e.g. Gummersbach 

et al. [19981 Rolleston et al. [2000t Daflon & Cunha[2004]) as 

Send offprint requests to: M.F. Nieva (nieva® sternwarte.uni- the basis for galactochemical evolution modelling (e.g. Hou et 

erlangen.de) al. 120001 Chiappini et al. 120011 Cescutti et al. l2006l l. Stellar 

* Based on observations collected at the European Southern Obser- data can hereby complement and independently verify analo- 

vatory, Chile (ESO N° 074.B-0455(A)) and at the Centra Astronomico gous abu ndance studies of H n regions (e.g. Shaver et al. [13551 . 

Hispano Aleman at Calar Alto, operated jointly by the Max-Planck Using ±Q present generation of large te lescopes and high- 
Institut fur Astronomie and the Instituto de Astrofisica de Andalucfa 



** Table [T]and Figs. U6H19l are only available in electronic form at 
http : //www . edpsciences . org 



resolution spectrographs, pristine abundances can be derived 
not only for Galactic early-type stars, but also for objects in the 
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metal-poor environments of the Magellanic Clouds (e.g. Korn 
et al. EU02ll2005l Rolleston et al. 120031 Hunter et al. [2Q06l >. 

In order to be meaningful, the parameter and abundance 
determination for individual objects has to be unbiased by sys- 
tematic error. Early and mid O-type stars pose considerable 
challenges to model atmosphere analyses, because of spheric- 
ity, mass-outflow, and non-LTE line-blocking/blanketing ef- 
fects. For these reasons the least luminous (but most numer- 
ous) massive stars, late O and early B-type (OB) stars (Jaschek 
& Jaschek 119901 ), have been at the focus of abundance stud- 
ies for a long time (e.g. Gies & Lambert 1992] Kilian[l992i 
Cunha & Lambert 119921 and numerous similar studies ever 
since). Overall, the atmospheres of OB dwarf and giant stars are 
supposed to be described reasonably well by one-dimensional, 
plane-parallel, homogeneous, and hydrostatic line-blanketed 
LTE models in radiative equilibrium. However, this does not 
imply that quantitative analyses of the spectra of OB dwarfs 
and giants are trivial. 

Most of the published quantitative analyses follow two ap- 
proaches. Either they are based on pure LTE analyses (e.g. 
Rolleston et al. (2000), and references therein) or they solve 
the restricted non-LTE problem on prescribed line-blanketed 
LTE model atmospheres. Typically, metal ions are considered 
in this case. A simultaneous treatment of H and He is excep- 
tional (e.g. Gummersbach et al. 119981 Morel et al. 120061 ). A 
few studies use unblanketed non-LTE atmospheres with subse- 
quent non-LTE line formation (e.g. Rolleston et al. 2000), and 
only recently have grids of fully line-blanketed non-LTE mod- 
els become available (Lanz & Hubeny 2003 concentrating on 
O-type stars but extending into the B star regime; Dufton et 
al. 120051 concentrating on B-type stars, not publically avail- 
able). A thorough test of the models regarding their ability to 
reproduce the H and He spectra via direct comparison with ob- 
servation throughout the entire visual range and in the near-IR 
in OB dwarfs and giants covering a broad parameter range has 
not been published so far (note their availability - at least for 
subsets of lines - for O stars: Bouret et al. 120031 Repolust et 
al. 120041 120051 Mokiem et al. 120531 [20061 Heap et al. 120061 
Having this kind of test available would help the user of grids 
to understand the strengths and the limitations of the models. 
Usually, only one or two hydrogen Balmer lines and selected 
helium lines in the optical blue are considered in the literature. 

Hydrogen and helium are of major interest in the astro- 
physical context, as they constitute practically all light-emitting 
plasma. The lines of hydrogen and helium are the strongest 
spectral features in OB stars. Inasmuch as they are primary di- 
agnostic tools for stellar analyses throughout the Hertzsprung- 
Russell diagram, they sample the plasma conditions through- 
out large parts of the stellar atmosphere, to a far greater extent 
than do the metal lines. However, the information on the tem- 
perature and density structure encoded in the spectra has to be 
interpreted. This is done by comparison with synthetic spectra, 
which requires that the basic atmospheric structure equations 
in combination with the radiative transfer problem be solved. 
The model predictions may differ, depending on the approxi- 
mations made and on the atomic data used. Their quality can 
be assessed by their ability to reproduce observation. In the op- 
timum case all observational constraints (continua/spectral en- 



ergy distribution, line profiles) should be reproduced simulta- 
neously, indicating the absence of systematic error (assuming a 
unique solution). A thorough reproduction of the hydrogen and 
helium line spectra should therefore be viewed as a precondi- 
tion for all further studies. 

Non-LTE effects play a dominant role in the formation 
of the hydrogen and helium line spectra in early-type stars, 
as known since the seminal work by Auer & Mihalas d 19721 
119731 ). Despite the enormous progress made over the past thirty 
years, some notorious problems have remained. Observations 
in the (near-)infrared provide one key to improving the sit- 
uation via extension of our observational database to a do- 
main of amplified non-LTE effects (in OB stars). Some of the 
problems have been related recently to the remaining inaccu- 
racies in the atomic data. Thus, the modelling of the hydrogen 
Paschen, Brackett, and Pfund series in early-type stars could 
be improved, resulting in corrections of equivalent widths by 
as much as a factor 2-3 (Przybilla & Butler 2004; but see 
also Repolust et al. 2005). Also, the observed behaviour of the 
He i /110830A transition in OB dwarfs could be reproduced 
for the first time (Przybilla 2005 ). In other cases, the reasons 
for shortcomings in the non-LTE modelling can be subtler. An 
example of this is the He i singlet line problem in early-type 
stars: computations with non-LTE model-atmosphere codes re- 
veal discrepancies not only between theory and observation but 
also between different theoretical calculations. The overlap of 
an He i resonance transition with Fe iv lines results in high sen- 
sitivity to the model assumptions (Najarro et al. 2006 ). 

The aim of this paper is to evaluate the status of non- 
LTE line-formation computations for the two most abundant 
elements in the most common targets of massive star analy- 
ses, OB dwarf and giant stars. This work constitutes the ba- 
sis for further studies of the metal line spectra (e.g. Nieva 
& Przybilla l2006bl [2007b . First, we test our hybrid non-LTE 
approach (Sect. 2) thoroughly on high-quality spectra of six 
stars in the solar vicinity (Sects. 3 & 4). In contrast to typ- 
ical studies from the literature, we investigate the entire hy- 
drogen and helium line spectra in the optical range, plus some 
additional near-IR data, taking advantage of our recently im- 
proved non-LTE model atoms. After making sure that excellent 
agreement between theory and observation can be obtained (i.e. 
also avoiding the aforementioned He i singlet line problem), we 
compare our modelling with libraries of synthetic fluxes from 
the literature (Sect. 5). This is done in order to test their suit- 
ability for quantitative analyses of OB dwarfs and giants. Such 
libraries are required for (automatised) analyses of large obser- 
vational datasets obtained with existing or future multi-object 
spectrographs (e.g. the VLT-FLAMES survey of massive stars: 
Evans et al. 120051 GAIA: Perryman et al. I2001I ). It is shown 
that reliable modelling of the line spectra of the two most basic 
elements is not straightforward. On the contrary, considerable 
systematic errors may result for quantitative analyses of OB 
stars when applying these libraries blindly. Finally, the results 
of our investigation are discussed (Sect. 6) and the main con- 
clusions summarised (Sect. 7). 
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2. Model calculations 

The hybrid non-LTE approach solves the restricted non-LTE 
problem on the basis of prescribed LTE model atmospheres. 
The approach is physically less elaborate than fully self- 
consistent non-LTE calculations, as more approximations are 
involved. However, at the same time it is superior to the pure 
LTE approximation. In particular, it provides an efficient way to 
compute realistic synthetic spectra in all cases where the atmo- 
spheric structure is close to LTE (which puts restrictions on the 
parameter space coverage). The hybrid non-LTE approach also 
allows extensive non-LTE model atoms to be implemented, fa- 
cilitating a highly detailed treatment of the atomic processes 
involved (e.g. account for the resolved resonance structure in 
photoionizations, avoidance of the - powerful, however also 
approximate - superlevel formalism). 

We compute line-blanketed, plane-parallel, homogeneous, 
and hydrostatic LTE model atmospheres using the ATLAS9 
code (Kurucz 1993b). Non-LTE population numbers and 
synthetic spectra are then obtained with recent versions 
of DET AIL an d SURFACE (Giddings HWTl Butler & 
Giddings I19851 I. The coupled radiative transfer and statis- 
tical equilibrium equations are solved with DETAIL, em- 
ploying the Accelerated Lambda Iteration (ALI) scheme of 
Rybicki & Hummer (1991). Synthetic spectra are calculated 
with SURFACE, using refined line-broadening theories. 

The non-LTE model atoms for hydrogen and He i/ii adopted 
in the present work are described in detail by Przybilla & Butler 
(2004) and Przybilla (2005), respectively. Use of improved 
atomic data for electron impact excitations, in particular from 
ab-initio computations, allows consistent results from the hy- 
drogen lines in the visual and near-IR to be derived through- 
out the entire range of early-A to O stars . A 15-level model 
is used for modelling main sequence stars, as well as a 20- 
level model for the giants. The Hei/n model atom considers 
all He i LS -coupled terms up to the principal quantum number 
n — 5 individually, and packed levels up to n = 8 separately 
for the singlet and triplet spin system. All levels up to n - 20 
are considered in the He n model. This model atom has been 
successfully used to reproduce observed trends of the highly 
non-LTE-sensitive He i A10 830 A transition in early-type main 
sequence stars (Przybilla 20051. The He i/n model was also em- 
ployed to analyse the visual/near-IR spectra of extreme helium 
stars (Przybilla et al. 2005) and subluminous B stars (Przybilla 
et al. 12006a* . 

Radiative transitions in DETAIL are treated with simplified 
line broadening formalisms: for transitions between hydrogen 
levels with n < 7 approximate Stark-broadening (Griem 1960 
following the implementation of Auer & Mihalas 119721 
Appendix) is considered, while for all other transitions, also 
in Hei/n, depth-dependent Doppler profiles are assumed. 
Microturbulence is explicitly accounted for by including the 
appropriate term in the Doppler width 

A^(^) 1/2 , (1) 

where A® is the rest wavelength of the transition, c the speed 
of light, o t h the thermal velocity for the chemical species of in- 
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Fig. 1. Sensitivity of theoretical Hei line profiles to modifi- 
cations of the microturbulent velocity in the statistical equi- 
librium calculations. The microturbulence of the ATLAS9 
model atmosphere structure is held fixed at ^ = 8kms _1 . The 
test calculations have been done for one of our sample stars 
(HR 3055). Similar effects are found for other He i lines, while 
the Hen and H lines are practically insensitive to this. 



terest and £ the microturbulent velocity. Both continuous ab- 
sorption and line blocking (via Kurucz' opacity distribution 
functions, ODFs, Kurucz fl993al using the 'little' wavelength 
interval versions) are accounted for as background opacities in 
solving the radiation transfer. In this regard the hybrid non-LTE 
approach has an advantage over present-day 'exact' non-LTE 
computations: all species responsible for metal line blocking 
and blanketing can be considered, though approximately. The 
'exact' non-LTE methods, on the other hand, are constrained to 
several abundant light elements and typically iron (and nickel), 
i.e. only the major line opacity sources are covered. Note that 
these ODFs were calculated for solar abundances according 
to Anders & Grevesse ( 119891 1. The latter have been revised in 
more recent work such as Grevesse & Sauval ( 1998). In partic- 
ular the reduction of the abundance of iron (the most important 
line opacity source) by ~0.2dex should be considered. This is 
done in our work by adopting the Kurucz dl993al ) ODFs for 
appropriately reduced metallicity; see Sect. l5.1l for further dis- 
cussion. 

The resulting non-LTE populations are then used to com- 
pute realistic line profiles with SURFACE. The same micro- 
turbulent velocity as in DETAIL and in the model structure 
computations with ATLAS9 is adopted. In this step of the cal- 
culation detailed Stark-broadening data are employed, as sum- 
marised in Table Q] (available in the online version). All other 
important data relevant to line formation are also given there: 
wavelengths, lower and upper levels involved in the transition, 
oscillator strengths log gf, their accuracies, and sources. 

Note that, in typical non-LTE computations for OB stars, 
microturbulence is only included in the final profile calcula- 
tion. Our choice is based on test calculations that indicate line- 
profile fits are improved if microturbulence is also included 
in computing the level populations. The net effect is a slight 
strengthening of the lines (Fig.[T). However, the effect is far less 
pronounced than described by McErlean et al. ( 1998), who in- 
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Table 2. Parameters of the programme stars 
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vestigated unblanketed non-LTE models for B-type supergiants 
at slightly higher microturbulence; cf. their Figs. 3 & 4. 

Our hybrid non-LTE approach involving the codes AT- 
LAS9, DETAIL, and SURFACE (henceforth abbreviated ADS) 
is tested here for early B-type dwarfs and giants, supplemented 
by LTE computations with ATLAS9 and SURFACE (AS). Our 
methodology may be applied to a wider range of stellar param- 
eters (i.e. effective temperature T e ff and surface gravity logg). 
Line-blanketed, static, and plane-parallel LTE models provide 
an even more realistic description of stellar atmospheres at 
lower temperatures and higher gravities (excited Hen states 
should be ignored at lower T e g in order to avoid numerical in- 
consistencies). Slightly higher temperatures (late O-types) and 
lower surface gravities (less-luminous supergiants) may also 
be covered, until the hybrid non-LTE approach meets its limita- 
tions when non-LTE effects on the atmospheric structure and/or 
hydrodynamic mass-outflow may no longer be neglected. 



3. Observational data 

We test our analysis technique on six bright Galactic objects in 
the entire optical range and also for near-IR lines when avail- 
able. The programme stars sample the parameter space in ef- 
fective temperature and surface gravity covered by typical ap- 
plications. 

High-S/N Echelle spectra of tSco (HR6165), HR3055, 
HR 1861, HR2928, HR3468, and HR5285 were obtained by 
M. Altmann using FEROS (Fiberfed Extended Range Optical 
Spectrograph, Kaufer et al. 1999) on the ESO 2.2m telescope in 
La Silla. The data reduction was performed within the FEROS 
context in the ESO MIDAS package, using optimum extrac- 
tion. The spectra were normalised by fitting a spline function 
to carefully selected continuum points. This suffices to retain 
the line profiles of the Balmer lines in these early-type stars. 
Finally, the spectra were brought to the wavelength rest frame 
by cross-correlation with an appropriate synthetic spectrum. Of 
the entire wavelength range covered by FEROS, only the part 
between ~3800 and 8000 A meets our quality criteria for fur- 
ther analysis. The spectra are compromised by the lower sen- 
sitivity of the instrument at shorter wavelengths, and the re- 
duced stellar fluxes in the far red. FEROS provides a resolving 
power R ^ A/AA « 48 000, with 2.2 pixels per AA resolution el- 
ement. An S/N of up to ~800 is achieved in B. With respect to 
resolution and signal-to-noise ratio, the spectra available to us 
are of much higher quality than in typical studies of OB stars, 
basically excluding observational uncertainties from the error 
budget. 
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Fig. 2. Scheme of the iterative procedure for obtaining a simul- 
taneous fit to the hydrogen and helium lines. 

A supplementary high-S/N spectrum of HR 1861, also cov- 
ering the higher Paschen series, was obtained using FOCES 
(Fibre Optics Cassegrain Echelle Spectrograph, Pfeiffer et 
al.[l998) on the Calar Alto 2.2m telescope. The data were pro- 
cessed in a standard way, using the data reduction routines 
described by Pfeiffer et al. ( [19981 . An £^40000 (2 pixels 
per AA resolution element) was achieved. In addition, a high- 
S/N spectrum in the Zf-band of r Sco taken with Subaru/IRCS 
(R a 12 000) is available for analysis; see Hanson et al. {2005) 
for details on the observations and data reduction. Finally, a 
high-S/N spectrum in the A 2.058 /urn region of r Sco taken with 
UKIRT/CGS4 (R 19000) is also available for analysis (Zaal 
et al. fT999l) . 

4. Applications to observations 

Theoretical profiles were fitted to observations in an iterative 
procedure summarised in Fig. 12 The final atmospheric param- 
eters T e g and logg, projected rotational velocities vsin/, and 
micro- and macroturbulent velocities (£, £) coincide with those 
derived in Nieva & Przybilla (2006a), which are further refined 
by use of the C ii/m and C ii/m/iv ionization equilibria (Nieva & 
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Fig. 3. Impact of stellar parameter variations on non-LTE pro- 
file fits, exemplarily for Hy, He 1/14026 A, and Heii/15411 A 
in HR 3055 (B0 III). Synthetic spectra for our final parameters 
(see Table [2] thick line) and varied parameters (thin lines for 
our uncertainty estimates and dashed lines for values typically 
found in the literature) are compared to observation. See the 
text for a discussion. 

Przybilla 2006b c), see Table [2] The impact of stellar parame- 
ter variations on non-LTE profile fits to Hy, He i /14026 A, and 
He n /1541 1 A in the hot giant HR 3055 is visualised in Fig. [3] 
Two values for the parameter variations are adopted, accord- 
ing to our uncertainties of 300 K/0.05 dex in r e ff/log g and typ- 
ical values from the more recent literature (lOOOK/O.lOdex). 
All other hydrogen Balmer and helium lines react in a similar 
way. The sensitivity of the hydrogen and He i lines to parameter 
variations is low, such that the uncertainties cannot be reduced 
much below the typical values even for high-S/N observations 
like ours. Only the He n lines are highly sensitive to changes in 
r e ff and, to a lesser degree, in log g. However, by taking metal 
ionization equilibria into consideration (e.g. C n/m or C ii/ra/iv), 
which are even more sensitive than He n lines, the parameters 
can be constrained more accurately. 

Projected rotational velocities, micro- and macroturbulence 
values have also been verified by fitting the carbon lines. 
Note that the comparatively high macroturbulence in HR 3055 
amounts to less than twice the sound speed in the atmo- 
spheric plasma. The macroturbulent broadening may thus be 
explained by a pattern of ascending and descending surface el- 
ements caused by (higher-order) nonradial oscillations (see e.g. 
Lucy 11 9761 1 . The microturbulence values are typically lower 
than found in previous work (e.g. Kilian 1992). Differences in 



J e ff and log g are also found. Solar helium abundances y (by 
number) are found in all cases. 

4.1. Visual 

Synthetic profiles for a selection of 6 hydrogen Balmer and 18 
He i/ii lines in the visual are compared with observation for the 
sample stars in Figs.|4]and|5]and in Figs.[T6HT9l(available in the 
online edition). They constitute our best simultaneous fits to the 
measurable H and He lines in the available spectra. Preference 
for our selection has been given to (mostly) unblended features 
with good broadening data. A summary of all available lines 
is given in Table Q] where blending species are also identified. 
Our hybrid non-LTE approach (ADS) allows us to reproduce 
the hydrogen Balmer and He i/ii lines in the visual more pre- 
cisely, with few (well-understood) exceptions. The ionization 
equilibrium of He i/ii puts tight constraints on T e ff in the two 
hotter stars. 

For the hottest dwarf of our sample, rSco, a very good 
match between model and observation is achieved (Fig.|4]l, ex- 
cept for the cores of Ho- and He n /14686 A. This is because of 
the neglect of the stellar wind; see Przybilla & Butler (2004) 
and Mokiem et al. (2005) for results of hydrodynamic com- 
putations. The discrepancies in Hei ,14121 A occur because of 
blends with metal lines (O n, C m, and Fe m; unaccounted for 
in our computations), which can be nicely resolved at this low 
v sin i. An improved fit to He i ,44921 A may be obtained with 
better broadening data for the forbidden component. The spec- 
tral region around Ha suffers from artifacts introduced by CCD 
defects that can only partially be compensated for in the data 
reduction process. 

Line fits to the hot giant HR 3055 are displayed in Fig. [5] 
Excellent agreement between theory and observation is also 
found in this case, with a significantly improved fit quality of 
Ha and Hen /I4686 A, because of an apparently weaker wind. 
This star shows a higher v sin i and ( than r Sco. Therefore the 
He i AA\2\ A blend is no longer resolved, leading to an appar- 
ently worse fit. 

He n /14686 A is the only visible (weak) feature of He n in 
the intermediate temperature stars HR 1861 and HR2928; see 
Figs.[T6]and[T7] Good fits are obtained for this line and the fea- 
tures of He i (i.e. establishing the ionization equilibrium) and 
hydrogen. Lines of He n are absent at even lower temperatures, 
HR 3468 and HR 5285 in our sample, see Figs. [Hand[T9] 

4.2. Near-IR 

Additional spectra are available in the near-IR for two stars. 
An excellent fit to the higher Paschen series is obtained for 
HR1861 in non-LTE, despite the relatively low S/N and the 
presence of telluric lines, see Fig. [6] Good agreement between 
the non-LTE spectrum synthesis and observation can also be 
obtained for the He i A 2.058 /im feature and practically perfect 
agreement for the A2.ll/um lines in tSco, see Figs. [7] and [8] 
Note that the A 2.058 fim transition is situated in an atmospheric 
window with a series of strong telluric lines, therefore its shape 
is highly sensitive to the detailed approach in the data reduc- 
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Fig. 4. Non-LTE line fits to observed hydrogen and helium features in t Sco (BO. 2 V), based on the atmospheric parameters 
summarised in Table [2] 
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Fig. 6. Modelling of the Paschen series of HR 1861 with our 
non-LTE (ADS) approach. Note the presence of numerous 
sharp telluric H2O lines. All theoretical spectra in the near- 
IR have been computed with the same atmospheric parameters 
than the models in the visual. 



Fig. 7. Modelling of the He 1 A2.058 jjm singlet line in r Sco. 



tion process (see e.g. Zaal et al. 1999). The two different obser- 
vations in Fig. [7] may exemplify the difficulty of accurate tel- 
luric line removal in the data reduction process (or alternatively 
an intrinsic time variability of the feature). Good agreement 
with our model is obtained in the case of the higher-resolution 
UKIRT/CGS4 spectrum. The LTE approach is not even able to 
reproduce the observation qualitatively. The line fits in the near- 
IR are based on the same atmospheric parameters (Table [2]) as 
used for the modelling of the optical spectra. 

Non-LTE effects can be amplified in the Rayleigh-Jeans 
part of the spectral energy distribution, as demonstrated in this 
case. See e.g. Przybilla & Butler (2004) for a discussion and 
for line fits to Bry in the /f-band and to additional Brackett and 
Pfund lines in this star. The case of Hei /U0830A has been 
discussed by Przybilla (2005). 




Fig. 8. Modelling of the Hei A2.ll /im singlet and triplet fea- 
ture in tSco. As also shown in Fig. [7] these near-IR transitions 
experience stronger non-LTE effects than the spectral lines in 
the visual. 



5. Comparison with other model predictions 

In this section we compare our hybrid non-LTE computations 
(ADS) with four other approaches. Two of them are available 
grids from the literature, and to understand their discrepancies 
to our ADS computations we calculate additional non-LTE and 
LTE models. 

I) We replace the ATLAS9 models by line-blanketed, 
plane-parallel, and hydrostatic non-LTE model atmospheres 
taken from the publically available OSTAR2002 grid (Lanz & 
Hubeny 2003 LH03) and carry out the non-LTE line formation 
calculation with DETAIL and SURFACE as described above. 

II) We use the non-LTE calculations from the OSTAR2002 
grid (TLUSTY and non-LTE line formation with SYNSPEC), 
published by LH03. 

III) We calculate LTE spectra based on ATLAS9 at- 
mospheres and subsequent LTE spectrum synthesis with 
SURFACE (AS). 

iv) We use the Padova grid (Munari et al. 120051 1. based 
on ATLAS9 atmospheres and LTE spectrum synthesis per- 



formed with the SYNTHE code (Kurucz & Avrett QMTJ 
Kuracz [T993cI 

5.1. Atmospheric structures, SEDs: LTE vs. non-LTE 

A comparison of LTE (ATLAS9) and non-LTE (TLUSTY) at- 
mospheric structures and spectral energy distributions (SEDs) 
computed with ATLAS9+DETAIL and TLUSTY is made in 
Fig. Models for a hot giant and a dwarf (T eB = 32 500 K, 
log g- 3.75 and 4.25, respectively) are considered, approxi- 
mately delineating the upper temperature boundary of our ob- 
servations (tSco). Reduced non-LTE effects can be expected 
for cooler models. 

Excellent agreement is found for the temperature and den- 
sity structures. This is a basic requirement for successful appli- 
cation of the hybrid non-LTE approach for spectrum synthesis. 
The temperature structures deviate by less than 1% in the in- 
ner atmosphere, including the regions where the weaker lines 
and the wings of the stronger features are formed (logm > -1 ; 
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Fig. 9. Upper panel: Comparison of ATLAS9 and TLUSTY temperature structures and electron densities (insets) as function of 
column mass. The computations have been performed for giant and dwarf models. Lower panel: Comparison of spectral energy 
distributions, the radiation field computed by DETAIL on the basis of the ATLAS9 atmospheric structure vs. TLUSTY. 
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Fig. 10. Comparison of the most discrepant hydrogen and 
Hei/n line profiles from our hybrid non-LTE approach (ADS) 
and a TLUSTY-DS calculation for a hot main-sequence model. 
Practically perfect agreement is obtained, with small discrep- 
ancies occurring only in the wings of He n /14686 A. 



see Sect. 15.3b . At the formation depths of the cores of the 



stronger H and He lines (-3 < log m < -1.5) the differences 
may increase to <2-3%. Stronger deviations may occur only 
in the outermost parts of the atmosphere, outside the line- 
formation depths. Note that this good a match is obtained 
only if the effects of metal line-blanketing are correctly ac- 
counted for. In particular, the computations should be made 
for identical metal abundances. This is complicated by the fact 
that the ODFs of Kurucz (1993$ were computed assuming 
scaled solar abundances from Anders & Grevesse ( 1989), while 
the TLUSTY computations assume abundances from Grevesse 
& Sauval dl9981 l. The most important difference is a down- 
ward revision of the iron abundance by ~0.2 dex in the later 
work. Consequently, we use ODFs with correspondingly 'sub- 
solar' metallicity in order to correct for the discrepancies in 
the line opacities, with [Fe/H] as a metallicity substitute. See 
also Przybilla et al. (2006b) for a discussion of such 'empiri- 
cal' corrections. We should note that while the differences are 
small at (near-)solar abundances the non-LTE effects on the at- 
mospheric structure will increase with decreasing metallicity. 
Nevertheless, our hybrid non-LTE methodology for OB star 
analyses should be applicable down to SMC metallicity, as in- 
dicated by an analogous comparison for models at 1/5 x solar 
abundances. 
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Fig. 11. Comparison of selected H and Hei/n line profiles from our hybrid non-LTE approach (ADS), non-LTE computations 
from TLUSTY+SYNSPEC, and two LTE calculations (ATLAS9+SURFACE and ATLAS9+SYNTHE) for a hot main-sequence 
model. 



The SEDs computed with ATLAS9+DETAIL and 
TLUSTY show excellent agreement over almost the entire 
wavelength range. Small differences occur in the EUV, most 
notably in the He n continuum. This is a significant improve- 
ment over the comparison of TLUSTY with ATLAS9 model 
fluxes (not shown here), which predict much lower fluxes 
in the Lyman and helium continua. The LTE computations 
neglect the non-LTE overionization of the hydrogen and He i 
ground states. This overionization reduces the bound-free 
continuum opacity; see Sect. 15.31 for a more comprehensive 
discussion. 



5.2. Spectra: hybrid non-LTE vs. full non-LTE and LTE 

Comparisons of synthetic profiles of several strategic lines of 
hydrogen and He i/ii are made for three test cases, where mod- 
els are available from the published grids. These frame the 
parameter space studied in the present paper. The test cases 
comprise a hot dwarf model (T e g = 35 000 K, log g- 4.5), at 



slightly higher T e tf than covered by our observation^ shown in 
Fig. QT] a hot giant model (32 500, 3.75), with similar temper- 
ature to tSco, see Fig. [TJl a cool giant model (20000, 3.00), 
with both r e ff and logg lower than covered by our observa- 
tions, Fig. [13] In order to be consistent with the published grids, 
our computations consider solar metal abundances (Grevesse & 
Sauval 1998) and solar helium abundance. Our synthetic spec- 
tra and those of the OSTAR2002 grid are degraded to the high- 
est resolution (R = 20 000) available from the Padova grid. Note 
that metal lines are neglected when we calculate the emergent 
spectrum, and they are considered indirectly via line blanket- 
ing/blocking effects. 

i) ADS vs. TLUSTY-DS. This comparison allows effects 
caused by differences in the model atmosphere structures to be 
disentangled. A practically perfect match for the (35 000,4.5) 
model (see Fig. [TDJ is obtained in the ADS and TLUSTY-DS 



1 We expect our approach to be valid at slightly hotter temperatures, 
at least in main sequence stars where the stellar wind does not influ- 
ence the photospheric layers strongly. 
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Fig. 12. Comparison of selected H and He i/n line profiles from our hybrid non-LTE approach (ADS), the non-LTE computa- 
tions with TLUSTY+SYNSPEC and our LTE approach ATLAS9+SURFACE for a hot giant model. Here, HR3055 may act as 
observational discriminator, indicating our results to be appropriate (see Fig.[5]l. 



computations, which share the same model atom. This indi- 
cates good agreement of the LTE and non-LTE atmospheric 
structures at even slightly higher temperatures than discussed 
in Fig. [9] We refrain from presenting further comparisons of 
these two approaches, as the discrepancies are even smaller at 
lower temperatures. 

n) ADS vs. TLUSTY+SYNSPEC (LH03). These results 
are obtained using two independent methods (model at- 
mospheres, model atoms, numerical solution). Nonetheless, 
good agreement is found on the whole for the (35 000,4.5) 
and (32500,3.75) models from an inspection of Figs. QT| 
and Q~2l respectively. Notable differences between ADS and 
TLUSTY+SYNSPEC occur in the line cores of the Balmer 
lines (the latter filled in by emission) and in the Hei sin- 
glet lines, which are systematically weaker in the case of 
TLUSTY+SYNSPEC, in contrast to observation (see Figs. H 
and |5]l. The Hei triplet lines derived from both approaches 
agree well. Small discrepancies occur in the line wings of the 



He 1 lines because of different broadening data. A good match is 
also obtained for the He 11 lines, with small differences arising 
in He H/14686 A. 

From the comparison of the TLUSTY+SYNSPEC and 
TLUSTY-DS results, which match our ADS, we can conclude 
that the aforementioned discrepancies arise because of subtle 
differences in solving the statistical equilibrium and radiative 
transfer problem. While our approach uses line opacities aver- 
aged over the ODF wavelength bins, the TLUSTY+SYNSPEC 
computations employ a more sophisticated opacity sampling 
technique. This however introduces a strong dependency on 
the model assumptions for Fe i\0 (a highly complex ion), which 
has lines overlapping with an He 1 resonance transition (Najarro 



2 In particular on the oscillator strengths of the Fe iv transitions in- 
volved. A reduction of the g/-values may alleviate the discrepancy 
between the He 1 singlet and triplet line strengths, seen, for example, 
in the OSTAR2002 grid. 
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Fig. 13. H and He i profiles for a cool giant model: our hybrid non-LTE approach (ADS) vs. our LTE (ATLAS9+SURFACE) and 
the corresponding Padova model (ATLAS9+SYNTHE); see Fig.[l8]for the closest observational analogue. 



et al. 2006 Jj. The same model atmosphere (TLUSTY) is used 
and both model atoms should be sufficiently robust for mod- 
elling the lines in the visual; see Przybilla & Butler (2004) and 
Przybilla (2005 ) for a discussion of this. 

m) ADS vs. ATLAS9+SURFACE. LTE computations with 
AS produce narrower Balmer lines for the (35 000,4.5) and 
(32500,3.75) models (the differences reducing progressively 
from Ha to the higher series members), which leads to over- 
estimated surface gravities in that case. At the same time, all 
He i lines are too shallow in LTE, the trend increasing from 
blue to red and showing larger discrepancies at lower grav- 
ity. On the other hand, rather good agreement is found for the 
He ii lines, the LTE predictions being slightly weaker than ADS 
for the hot giant. In Fig. [13] a comparison of our hybrid non- 
LTE with our pure LTE prediction is made for a (20000,3.0) 
model with a temperature slightly below than the lower limit 
of our programme stars, and at significantly reduced surface 



3 As a further test we have calculated synthetic spectra for 
the (35 000,4.5) and (32 500,3.75) models on the basis of unified, 
line-blanketed non-LTE model atmospheres (FASTWIND, Puis et 
al. 1200 5 ). Excellent agreement with our ADS results is found for both 
the He i singlet and triplet lines. Note that FASTWIND also uses an 
approximate treatment of line blocking. 



gravity. Here, the wings of the Balmer lines show much bet- 
ter agreement than at higher temperatures (cf. Fig. ITTb. as well 
as the He i AA 4437 and 4713 A lines. The line cores are also 
discrepant, increasingly so from H6 to Ha. Many of the He i 
lines experience significant non-LTE strengthening, in partic- 
ular those in the red. The line broadening data is the same in 
ADS and AS, so the He i wings are very similar. The forbidden 
components are also accounted for in both approaches. 



iv) ADS vs. ATLAS9+SYNTHE (Munari et al. I2005I 
The differences of these approaches were quantified for the 
(35 000,4.5) model, when possible. The Balmer lines from the 
Padova model present similar characteristics as the LTE AS ap- 
proach (in Fig.QT]they coincide), resulting in lower equivalent 
widths by up to ~30% relative to ADS. When using the Hy 
wings as a surface-gravity indicator, this translates to a system- 
atic error in log g by -0.2 dex, implying even larger errors for 
fits to the Hf3 and Ha wings. The Hei lines are generally too 
weak, by up to a factor of more than 2 in equivalent width, 
and the He n lines too narrow. For the most part, these discrep- 
ancies stem from the neglect of non-LTE effects on the line- 
formation process, as the differences in the atmospheric struc- 
tures are practically insignificant. 
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Fig. 14. Departure coefficients fo, of some strategic hydrogen and helium levels as a function of Rosseland optical depth Tr oss . 
The comparison is made for the giant (left) and dwarf atmospheric models (right column) already discussed in Fig. [9] for our 
approach (ATLAS9+DETAIL, thick lines) and the results of LH03 (TLUSTY, thin lines). Each level is coded by different line 
styles; see the legend in the corresponding panels. Line-formation regions (from core to wing) corresponding to our calculations 
are indicated. See the text for further discussion. 



Another limiting factor of the ATLAS9+S YNTHE compu- 
tations is the use of insufficient Stark broadening data (Voigt 
profile with constant Stark damping parameter). Our AS and 
ADS approaches improve on this, as realistic broadening data 
is used (see Table Q). In the ATLAS9+SYNTHE approach, it 
will not be possible to obtain reasonable agreement for the He i 
and Hen spectra at the same time. For the (20000,3.0) model, 
the Padova profiles are more similar to the AS approach. 
However, the diffuse He i lines still suffer from inappropriate 
broadening data, in particular the forbidden components are 
unaccounted for. The Hei lines are affected by non-LTE 
strengthening, increasing to the red. Only few Hei lines are 
quite similar in the three approaches: AA4437 and 4713 A 
match quite well, as do /L13867, 4121 (despite blends with 
metallic lines) and 5015 and 5047 A, not displayed here. 

We should mention that the published libraries of synthetic 
spectra were computed with different values of microturbulent 
velocity (OSTAR2002: 10km s 1 ; Padova: 2 km s 1 ). The ADS 



and TLUSTY-DS calculations with £ = 10 km s _1 were adopted 
for the comparison in Fig. QT| We made tests with a reduced 
£ = 2kms~ 1 , resulting in only small changes in the Hen pro- 
files - the most sensitive to modifications of The differences 
between the Padova grid and the other approaches are indeed 
due to the neglect of non-LTE effects and to additionally insuf- 
ficient broadening data and not because of discrepant microtur- 
bulences. The ADS and AS computations in Figs. [T2land[T3l 
were performed using the £ of the respective libraries. 

5.3. Line formation: hybrid non-LTE vs. full non-LTE 

We now try to identify the physical reasons for the differences 
in the non-LTE line profiles of hydrogen and helium in the last 
comparisons by a closer study of the underlying line-formation 
processes. For this we investigate non-LTE departure coeffi- 
cients and line source functions for three representative hydro- 
gen and six Hei/n lines, as derived in our ATLAS9+DETAIL 
and the TLUSTY computations. We choose the same models 
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Fig. 15. Ratio of line source function Si to Planck function B v at line centre as a function of tr oss for selected spectral lines of 
hydrogen and helium. The comparison is made in analogy to Fig. [14] for our approach (ATLAS9+DETAIL, thick lines) and 
LH03 (TLUSTY, thin lines). The individual spectral lines are encoded by the different line styles indicating the line-formation 
depths (from our calculations). See the text for further discussion. 



as discussed in Fig. [9] For the (32500, 3.75) model, a direct 
comparison with the resulting line profiles is facilitated by in- 
spection of Fig. [T2] 

For the levels involved in the transitions of interest and the 
hydrogen and helium ground states, departure coefficients £>, 
(referred to the ground state of the next higher ion) are dis- 
played in Fig. [14] The non-LTE effects on the level occupations 
give rise to departures of the line source function S\ from the 
Planck function B v ; see Fig.[l5]for a comparison of S\/B v from 
the ATLAS9+DETAIL and TLUSTY computations. We recall 
that 



bi/bje\p(hv/kT) - 1 

with h the Planck and k the Boltzmann constant, v the transi- 
tion frequency, and T the local temperature. For given T e g the 
non-LTE effects are strengthened with decreasing surface grav- 
ity, implying lower particle densities and thus larger mean-free- 
paths between photon absorptions. 



Hydrogen. We study three hydrogen lines, Ha, HB and Hy. 
Departure coefficients for levels n > 5 behave similarly to those 
for n = 5, which already traces the behaviour of the continuum 
closely. Consequently, the line source functions for the higher 
Balmer lines are similar to that of Hy. The non-LTE depopula- 
tion of the H ground state (Fig. [141 reduces the Lyman contin- 
uum opacity, giving rise to higher EUV fluxes than in LTE. The 
Lyman lines are expected to experience non-LTE weakening. 
Note that the TLUSTY calculations indicate a slightly stronger 
non-LTE depopulation of the ground state at continuum forma- 
tion depths than in our case. Good agreement of the departure 
coefficients for n = 2 is found, which is overpopulated at line 
formation depths. In combination with the higher H states be- 
ing close to LTE, this explains the non-LTE strengthening of 
the Balmer lines. At the formation depths of the line cores, in 
particular for Ha, the TLUSTY results show a less pronounced 
overpopulation, eventually leading to an underpopulation of the 
n = 2 state in the outer atmosphere. This explains the shallower 
lines from the TLUSTY computations relative to ADS, which 
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is a consequence of the upturn of S\/B Y (Fig. IT3T>. The effect 
decreases from Ha to the higher Balmer lines, as the core for- 
mation depths shift to deeper atmospheric layers. 

He i Singlets. We investigate two representative features: 
/L14921 and 6678 A. In general, the departure coefficients for 
most of the levels show rather good agreement, in qualitative 
behaviour as well as quantitatively. Notable differences in the 
ground-state overionization occur in the outer atmosphere. 
More relevant is the behaviour of the 2p 1 P° level, the lower 
level of practically all Hei singlet transitions in the visual. 
The non-LTE overpopulation in the TLUSTY results is far 
less-pronounced than in our case (Fig.[T4b. in particular for the 
giant model. This gives rise to discrepant line source functions 
(Fig. rT3T > and consequently differing line profiles (Fig. [T2] ) 
in the two approaches, with the ADS computations correctly 
predicting the non-LTE strengthening and thus reproducing 
observation (Sect. 0}. The states at higher excitation energies 
(«>4) are in LTE relative to the Hen ground state at line- 
formation depths. 

He i Triplets. We investigate two representative lines: AA4471 
and 5875 A. The departure coefficients from the two ap- 
proaches differ only slightly at line-formation depths (Fig.[T4"li. 
As a consequence, the source functions (Fig. [13] ) are also 
similar, resulting in negligible differences of the line profiles 
(Fig. fT2l >. Again, the states with n> 4 are in detailed balance 
relative to the He n ground state at depths relevant for the line 
formation. 

Hen. Two features were studied: /U4541 and 4686 A. The 
higher He n levels are close to LTE with the continuum state 
at line-formation depths. Only the n — 3 level shows a rele- 
vant overpopulation, resulting in a non-LTE strengthening of 
the ^4686 A line. The differences in the ADS and TLUSTY 
departure coefficients give slightly shallower profiles for this 
line in the OSTAR2002 model (Fig. [T2l . The source function 
for A 4541 A is essentially Planckian in the relevant region. 

6. Discussion 

Quantitative spectroscopy has to avoid systematic errors in 
order to be meaningful. A correct interpretation requires si- 
multaneous reproduction of all observational constraints (con- 
tinua/spectral energy distribution, line profiles). The line spec- 
tra of the most abundant elements, hydrogen and helium, are 
particularly important in this. 

We have investigated the suitability of hybrid non-LTE line- 
formation computations for quantitative analyses of the hy- 
drogen and helium line spectra of OB dwarf and giant stars. 
Our computations simultaneously reproduce the line spec- 
tra throughout the visual and near-IR (where available) at 
high quality, as well as the measured spectral energy distri- 
butions from the UV to near-IR, as discussed by Nieva & 
Przybilla (2006a). The only exceptions in our observational 
sample are the cores of Ha and He n /14686 A in t Sco, be- 
cause our computations do not account for stellar winds. For 



two He i lines blueward of the traditionally analysed spectral 
region (~4000-5000A), appropriate Stark broadening data is 
unavailable at present (see TableQ]). 

Comparisons with other approaches, including the repre- 
sentative libraries of synthetic spectra computed in non-LTE 
(OSTAR2002 grid of Lanz & HubenyUm and LTE (Padova 
grid by Munari et al. 120051 1 were also made. We find good over- 
all agreement between hybrid and full non-LTE calculations 
within the parameter space investigated here: for the model at- 
mosphere structures, spectral energy distributions, and many 
detailed line profiles. Our hybrid approach is able to simulta- 
neously reproduce the He i singlets and triplets, confirmed by 
comparison with observations. This is not trivial and becomes 
important in particular when only few He i lines are available 
in the observed spectra. 

A comparison of TLUSTY+SYNSPEC with other state- 
of-the-art non-LTE codes (F ASTWI ND and CMFGEN) for O- 
type stars by Najarro et al. (2006) implies differences in all 
modelling results for the Hei singlets. However, while moti- 
vated by the problems with accurately reproducing observa- 
tions, this study does notpresent a detailed comparison of mod- 
els and observed spectraj. Our work indicates that the He i sin- 
glet problem also persists in the hotter early B-type stars at 
solar metallicity (note that a significant reduction of metallicity 
alleviates the problem, as Heap et al. (2006) find a simultane- 
ous match of singlets and triplets in SMC O stars). 

Computations in LTE from the Padova grid, on the other 
hand, systematically predict too shallow and/or too narrow line 
profiles. In particular, the differences in the Hy wings - a com- 
mon surface gravity indicator - result in systematically overes- 
timated gravities by up to ~0.2 dex in LTE (for fixed T e ff ). The 
differences in the equivalent widths of the H lines can amount 
to up to ~30% and in the He i/n lines up to a factor >2 compared 
to our non-LTE calculations, with the discrepancies increas- 
ing with effective temperature. Nevertheless we are not able to 
quantify the differences in effective temperature determinations 
from non-LTE and LTE ionization equilibria of He i/ii, as some 
of the profiles of the Padova grid do not reproduce observations 
even qualitatively. 

In terms of parameter range and the underlying physics, 
the hybrid non-LTE approach is certainly restricted. It may 
be of limited use at higher temperatures (early and mid-O- 
type stars), lower gravities (early B-type and O-type super- 
giants), stars with strong winds, or extremely low metallici- 
ties. Nevertheless, the hybrid non-LTE approach is sufficient 
for studying normal OB dwarfs and giants, as it allows the ob- 
served line spectra to be reproduced in the visual and near-IR 
over a wide range of atmospheric parameters. Here it has ad- 
vantages over other more sophisticated non-LTE techniques: 
I) allows highly robust and detailed model atoms to be imple- 
mented and to be tested efficiently (i.e. concentration on atomic 
data while avoiding further complications like stellar winds), 
e.g. for metals with hundreds of levels and thousands of transi- 
tions; and II) the model calculations are fast: the computation 



4 Our own tests indicate that both the ADS calculations and 
FASTWIND provide high-quality fits to observation of the Hei sin- 
glet and triplet spectra in the parameter space of interest here. 
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of one H & He i/n model with DETAIL + SURFACE takes only 
a few minutes on a modern PC (as of 2006). 

7. Conclusions 

A hybrid non-LTE approach, such as is often employed in the 
literature for analyses of metal line spectra, has been thor- 
oughly tested for hydrogen and helium line formation in OB 
stars for the first time, using recently improved model atoms. 
The synthetic spectra simultaneously match almost all measur- 
able hydrogen and helium lines in the optical and (where avail- 
able) also in the near-IR spectral range of six test stars covering 
a wide range of stellar parameters. 

The comparison of state-of-the-art line-blanketed non-LTE 
and LTE models confirms that the atmospheric structure of 
OB dwarf and giant stars is described well under the as- 
sumption of LTE, but not their spectral energy distribution 
and also not their line spectra. For these stars in the range 
20000K<r eff <35 000K and 3.0 < log g< 4.5 (far from the 
Eddington limit), our hybrid non-LTE approach is equivalent 
to full hydrostatic non-LTE computations. 

It succeeds also in providing synthetic spectra that correctly 
reproduce the observed He i singlet lines, avoiding inconsisten- 
cies recently reported in the literature. 

In contrast to this, pure LTE modelling based on the Padova 
grid (or equivalent computations with ATLAS9+SYNTHE) 
may give rise to considerable systematic errors in the atmo- 
spheric parameter determination (T e ff, log g) and to subsequent 
elemental abundance studies for the hotter stars in particular. 
The problems with this approach may be remedied at least near 
the lower temperature boundary by implementing proper line 
broadening data. 

Finally, we suggest that the Hei triplet be preferred for 
analyses over the singlet lines from the OSTAR2002 grid in 
order to avoid systematic uncertainties, in agreement with the 
findings of Najarro et al. (2006). Hydrogen and helium line 
profiles from the Padova grid may be useful for quantitative 
analyses of stars cooler than ~22 000 K. Preference should be 
given to the He i ^3867, 4121, 4437, 4713, 5016, and 5048 A 
transitions, which are least affected by non-LTE effects. 
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Table 1. Atomic data for H and He i/ii line formation in the visual and NIR 
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a unpublished, priv. comm.; b vacuum wavelengths 



accuracy indicators - uncertainties within: AA: 1%; A: 3% 

sources of g/-values: FTS: Fernley et al. {I987l >; GRC: Green et al. {1957}; SPA: Schiff et al. {I97~fl 

sources for Stark broadening parameters: BCS69: Barnard et al. d!969l >; C: Cowley ( 11971b ; DSB: Dimitrijevic & Sahal-Brechot ( 1990); G60: 
Griem < fl960l ; G64: Griem {19641 ; G67: Griem {19671 ; GBKO: Griem et al. {19621 ; S: Shamev {19691 ; SB: Schoning & Butler {19891 ; SH: 
Stehle & Hutcheon {1999) 
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Fig. 16. Line fits for HR 1861 (Bl IV). For atmospheric parameters^sei=Airable[2l and for further discussion see the teWfi 1 86 1 
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Fig. 17. As Fig. Q2] but for HR2928 (Bl IV). 
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18. As Fig.dl but for HR3468 (B 1.5 III). 



HR 3468 
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. 19. As Fig.HU but for HR5285 (B2 V). 
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